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Abstract: In this publication, an algorithm is presented that combines the ME+PS approach 
to merge sequences of tree-level matrix elements into inclusive event samples [I] 
with the POWHEG method, which combines exact next-to-leading order matrix 
element results with the parton shower [3J. It was developed in parallel to the 
MENLOPS technique discussed in [3] and has been implemented in the event 
generator SHERPA [4] . The benefits of this approach are exemplified by some first 
predictions for a number of processes, namely the production of jets in e + e~- 
annihilation, in deep-inelastic ep scattering, in association with single W, Z or 
Higgs bosons, and with vector boson pairs at hadron colliders. 



1 Introduction 

In the past decade, the incorporation of higher-order corrections into parton-shower simulations has been in 
the centre of formal improvements of existing general-purpose Monte-Carlo event generators like HERWIG [5], 
PYTHIA [B], and SHERPA [3]. As parton showers approximate higher-order matrix elements only in the soft 
and collinear limits of the real-emission phase space, they do inherently not yield a good prediction in 
the hard domain. This deficiency can lead to rather large discrepancies between Monte-Carlo predictions 
and experimental data and should therefore be corrected. In this endeavour, two somewhat orthogonal 
approaches have been pursued. 

To improve the description of hard QCD radiation, "merging algorithms" (ME+PS), have been proposed [TJ, 
which were shown to be correct up to next-to-leading logarithmic accuracy in e + e~-annihilation into hadrons. 
These methods were improved and extended in a series of publications 8, 9 . A reformulation which generi- 
cally maintains the logarithmic accuracy of the parton shower was achieved in [TQ providing also the means 
to classify the various other methods and implementations according to their formal accuracy. It makes use 
of "truncated" parton showers, which are a concept initially introduced in [5] in the context of the POWHEG 
method. An alternative merging technique, known as "MLM merging" , was suggested in |12j and described 
in more detail and compared with the other algorithms in [13 . 

While the ME+PS methods succeed at improving the simulation of multijet events, they do not address the 
apparent problem that the total cross sections of the simulated inclusive samples are still of leading-order 

1 The algorithm has been further extended to include QED effects 1101 and multi-scale problems lll| . 



accuracy. However, this accuracy often is insufficient for tests of physics within and searches beyond the 
Standard Model. Examples for such situations range from luminosity measurements at the LHC through 
the production of W bosons to the determination of Yukawa couplings of the Higgs boson, once it is found. 
There has been a number of proposals of how to include at least the full next-to-leading order results in 
the simulation, but only two of them have been fully worked out and implemented in publicly available 
programs. 

The first method, Mc@NLO, pioneered in [14], and implemented for a large number of processes [15] relies 
on using the parton-shower kernels and their universal behaviour in the soft and collinear limits to sub- 
tract the infrared divergences of real-emission contributions to the NLO cross section. The parton shower 
then starts from either a Born-like configuration or from a configuration determined by the residual real 
correction contribution of the NLO calculation. By construction, there is some dependence on the details 
of the actual shower algorithm, which, to a certain extent, up to now seemed to limit the versatility of the 
method. This dependence was overcome by the second method, POWHEG, which was initially presented 
in [2J. This technique essentially is an improvement of an ingenious reweighting method, known for nearly 
two decades and applied individually to a plethora of processes |16j . To promote it to full NLO accuracy, 
it was supplemented with a local, phase-space dependent Zf-factor. The POWHEG method has been worked 
out for a number of processes |17j . using different parton-shower algorithms. A framework incorporating the 
core technology, independent of the specific parton-shower implementation and the matrix elements for the 
processes in question has been published in [TB]. It should be noted, though, that the POWHEG method relies 
on truncated parton showers in order to maintain the logarithmic accuracy in the radiation pattern, and it 
is this constraint, which introduces a parton-shower dependence. This has been an open issue and it was 
discussed in some of the above works, since older Fortran versions of PYTHIA and HERWIG typically do not 
provide truncated parton showers. Although it was typically assumed that the omission of truncation does 
not lead to sizable effects, the findings in [9] suggest that this depends significantly on the parton-shower 
model and the respective mismatch of its ordering parameter with the definition of the hardness scale. 

Having at hand two, somewhat orthogonal, methods (ME+PS and POWHEG) to improve both the hard QCD 
radiation activity and the total event rate in a given process, the question naturally arises whether it is pos- 
sible to combine both into an even more powerful approach. This is the topic addressed by this publication, 
resulting in a practical algorithm for merging both techniques. This algorithm has been implemented into 
the multi-purpose event generator SHERPA [3] , and will be made publicly available in a future version of the 
program. In a parallel development, Hamilton and Nason [3] suggested an identical method; however, their 
actual implementation only approximates the formal result. Due to the formal equivalence of both proposals, 
we will generally refer to the new technique as the MENLOPS approach. For a detailed presentation of the 
method, the reader is referred to Sec. [3] while Sec.[2jis used to introduce the notation employed in this paper. 
In Sec. [4] some first predictions of the SHERPA implementation are exhibited, exemplifying the improvements 
that can be achieved. These results and their interpretation solidify the findings of |3J. Section [5] contains 
our conclusions. 

2 Improving parton showers with higher-order matrix elements 

This section is devoted to the introduction of a formalism that allows to describe, on the same footing, the 
two basic methods for correcting parton-shower algorithms with real-emission matrix elements. A common 
notation will be established, which will later allow to merge the POWHEG method as reformulated in [TDJ, with 
the ME+PS method proposed in [T] . Like in the discussion is restricted to a dipole-like factorisation [20] . 
While the ME+PS approach allows to incorporate matrix-element information for arbitrary final states, in 
the POWHEG method only the first emission is fully corrected. To compare, and, ultimately, to combine both 
algorithms, only the expressions for the differential cross section describing the first emission off a given core 
interaction must therefore be worked out. This is where the combination takes place. 

Denoting sets of n particles in a 2 — > (n — 2) process by {a} — {a\, . . . , a n }, while their respective flavours 
and momenta are specified separately through {/} = {/i, ■••,/«} and {p\ = {pi,...,p n }, the generic 
expression for a fully differential Born-level cross section can be written as a sum over all contributing 
flavour combinations: 




where 




(2.1) 



{/} 
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The individual terms in the sum are given by 

B({a}) = £({a})B({a}) , B({a}) = — ^ — * \M B \ 2 ({a}) , 

(22) 

d<MM) = ^i^d^(M) , £(M;a* 2 ) -^i/ /i (x 1 ^ 2 )x 2 / /2 (x 2i/ i 2 ) , 

X\ X2 

where |A^b| 2 ({a}) denotes the partonic matrix element squared, d<?s({a}) is the corresponding differential 
n-particle partonic phase-space element, <?({/}) is the symmetry factor, F({a}) is the flux factor, and C is 
the parton luminosity. In the case of leptonic initial states, and ignoring QED initial-state radiation, the 
parton distribution functions f(x,fi 2 ) are replaced by <5(1 — x). 

In a similar fashion, the real-emission part of the QCD next-to-leading order cross section can be written 
as a sum, depending on parton configurations {a\, . . . , a„+i}, by replacing the Born-level matrix elements 
B with the real-emission matrix elements 1Z and the Born-level phase space d$s with the real-emission 
phase-space d&n. 

Furthermore, it is useful to introduce a notation for mappings from real-emission parton configurations to 
Born- level parton configurations and vice versa. They are given by (cf. |19j ) 

{p}^{0} and r*,kU«**\B,i*})-\ {£}^{p} ■ (2 - 3) 

The map bij k^a) combines the partons en and a,j into a common "mother" parton a^, in the presence of 
the spectator a& by defining a new flavour and by redefining the particle momenta. The inverse map, 
r~ fc(a) determines the parton configuration of a real-emission subprocess from a Born parton configuration 

and a related branching process ij, k — > ij, k. The radiative variables &r\b are thereby employed to turn the 
n-parton momentum configuration into an (n + l)-parton momentum configuration. 

Following [19], the real-emission matrix elements, lZ({a}), can be decomposed into a number of terms ~R-ij,k 
as 

n ihk {{a}) := PtJ A{S})K({a}) , where ftj ,,({a}) = ^ ^vAW}) (2 4) 

where are the dipole terms defined ibidem. Therefore, the real-emission differential cross section can 

be rewritten as a sum of trivially factorised contributions 

da R ({S}) = E E da B (b^ k ({d}))da l ^ B ({a}) , (2.5) 

{i,j} kjti,j 

where 

and where I Wfc = £7£y-,fc, cf. [19]. The no-branching probability of a parton-shower algorithm is derived, 
based on Eq. (2.6 1 and the decomposition of the dipole terms 2?jj & in the parton-shower approximation. It 
reads 

A< PS ) (*',*";{<?}) = J] AgJ (*'>*";{«})- ( 2 -7) 

where the dipole-dependent no-branching probabilities are given bj|^] 

j-''dl '"■ ■ rl:: 

fi=q,g 



(2. 



1 a s £{r~~ k {fi,t,z,4>;{a})]t) 

x s^*^''*'^ 



2 Note that, if initial-state partons are involved in the splitting, aL P ?' is a constrained no-branching probability, as x\ in < 1 
before and after the branching process. 
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The quantities Kij,k are the parton-shower evolution kernels, which depend on the parton flavours ft, fj and 
fk and on the radiative phase space ^^j^j = {t, z, </>}, cf. |19j . The variable t, with t oc 2piPj, is the evolution 
parameter while z is called the splitting variable of the parton-shower model and Jy^ is the Jacobian factor 
associated with the transformation of integration variables. 



2.1 The POWHEG approach 



Using a simple corrective weight, it is possible to modify the parton-shower such that it produces the 0(a s ) 
radiation pattern of the matrix element [19J. The corresponding dipole-dependent no-branching probability 
reads 



AgV,f";{a}) =ex P . 



E 

fi=Q,9 



16tt 2 



t" 



dt 



dz 



— J i:j<h (t,z,4>) 



1 S(r 5 ,s(/<;{/») ^Ar^ k {h,t,z^;{a})) 



(2.9) 



Sij S({f}) B({a}) j ' 

where the actual matrix elements R and B, which also include the respective parton luminosity factors, 



replace the parton shower kernel K. and the parton luminosities of the equation above, Eq. (2.8| 



The key point of the POWHEG method is, to supplement Monte-Carlo event samples from such matrix-element 
corrected parton showers with an approximate next-to-leading order weight to arrive at NLO accuracy. This 



is achieved by multiplying Eq. (2.9) with the local if -factor B({a})/B({a}), where 
B({a}) =B({a})+V({a})+I({a}) 

%,fc( r rj,fc({«})) - s y,fc( r j5,fe({«})) 



E E 

{ij.k} fi=i,g 



(2.10) 



In this expression, V({a}) is the NLO virtual contribution, including the collinear mass factorisation coun- 
terterms, while Sy,fe(r~ %(a)) and I({a}) denote real and integrated subtraction terms, respectively. 
This yields the following master formula for the value of an infrared and collinear safe observable, O 



(O) 



(POWHEG) 



E / d<MM)B(M) 
{/} 



A( ME )(t , M 2 ;{a}) 0({a}) 



unresolved 

fl 2 /"Zmax /'27T 



dt 



dz 



E E 16tt 2 /, 

{ij,k} h=q,g 

1 S(r r ~ k ({f})) Ry,fe(r r7 i({a})) 



— Jij, k (t, z, (f) 0(r~ ^({a})) 



(2.11) 



S « S({f}) 



B({a}) 



resolved 



Note that the second term in the square bracket describes resolved emissions, simulated by the matrix-element 
corrected parton shower, while the first term incorporates unresolved emissions and virtual corrections. 



To reveal the fixed order properties of (2.11 1 it is useful to inspect its expansion keeping terms up to 0(a s 



{0) (powheg) = gy d<MM) (b+V + i) ({S})0({a}) 

R«,fe(^,£({a})) e(*o- *) - S iJlfc (r- S ({S})) 



{/} 



0({a}) 



EE E /d*B(W)d*g* 

{/} {Q,fc} /i=?.9 

E E E / d$ B(M)d$^%, fc (r r ^({^))e(t-to)0(r^({a})) + G(a 2 ) . 



(2.12) 



{/} {«>fe} /*=9'9 
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It is imperative to note that neither the presence nor the precise form of A^ ME - 1 in the resolved emission 
term influence the fixed-order accuracy of the method at 0(a s ). Similarily, different choices of scales at 
which a s is evaluated in R, the B-function, and the matrix-element corrected parton shower emission terms 
contribute at 0{a 2 s ) because a s (/ii) = a s (/i2) (1 + 0(a s )) [3J. 

In order to obtain the correct leading logarithmic behaviour of the real-emission term, a s and the parton 



luminosities in A^ ME ) and the resolved emission term of Eq. (2.111 must be evaluated at scale k\. NLL accu 



racy can be restored for processes with no more than three coloured partons by means of the replacement |21j 



1 



2tt 



07 
18 



7T 

IT 



C A 



9 



(2.13) 



where the MS expression of a s should be used. In our Monte-Carlo implementation of the POWHEG method, 
we follow this approach. 



2.2 The ME+PS approach 

The ME+PS approach for the inclusion of matrix-element corrections into the parton-shower relies on a 
twofold generation of radiative corrections: Through an ordinary parton shower on the one hand and through 
real-radiation tree-level matrix elements on the other hand. In contrast to the POWHEG method, which only 
corrects the first emission off the core interaction, the ME+PS technique can be employed for arbitrary 
higher-order tree-level configurations. 

A first algorithm to achieve this was presented in [7j . The solution there is based on separating the radiative 
phase space into a region of soft /collinear emissions, the parton-shower (PS) region, and a region of hard 
emissions, the matrix-element (ME) region. By demanding each region to be filled by the respective way of 
generating radiation and some reweighting double counting and other problems can be avoided. The original 
formulation has been tremendously improved in [T], by realising that, formally, separate splitting kernels in 
the ME and PS regions can be defined, which add up to the full splitting kernel: 



t-ME 

ij,k 



(t, z, 4>) = )C ijtk (t, 2, 4>) e [Qij^it, z, 4>) - Q cut ) 
= JCij, k (t, z, 4>) e (Q cu t - Qi 3 ,k{t, z, 4)fj 



(2.14) 
(2.15) 



The functional form of the separation criterion Qij^ is in principle arbitrary as long as it identifies soft and 
collinear divergences in the real-radiation matrix elements. The approach, fully outlined in pQ, then replaces 
the splitting kernels in the ME region by the ratio of the real-emission and Born matrix elements, just like 
this is done in a matrix-clement corrected parton shower. However, in contrast to a reweighting technique, 
only emission terms are modified and no correction is applied to the no-emission probabilities. 

The ME+PS technique can be implemented in a master formula for the first emission, describing the expec- 
tation value of an arbitrary infrared safe observable O, similar to the POWHEG case: 



(O) 



(MEPS) 



J d<MM)B({a}) 

{/} 



A( PS )(i , M 2 ;{a}) 0({a}) 



unresolved 



E E 16^2 

{i],k} fi=q,g 



At 



dz 



2tt 



Jij, k {t,z,4,) 0(r~ k ({a})) 



1 8ira s C(r rk (f t ,t,z,<j>;{a});t) 

X — I : fcij,k{t,Z,(j)) , 6 I Qcut - Qij,k{t, Z, 



S, 



t 



C({a};t) 



(2.16) 



resolved, PS domain 

%r,,fe({/})) RiJ,*(r«is({3})) 



S({f}) 



B(M) 



e 



resolved, ME domain 



There are three components to the differential cross section: The term describing unresolved emissions, which 
is generated in the standard parton-shower approach, and the resolved part, which is now split between the PS 
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and the ME domain. Within the ME domain, the matrix-element generator is directly invoked to define the 
real-emission configuration. This is possible due to the restricted phase space, removing all infrared divergent 
regions by applying the cut in Qij^ and rendering the matrix clement finite. In this case, the Sudakov form 
factor A( PS \ which makes the matrix element exclusive, must be added explicitly. It can either be calculated 
analytically, like in the original formulation of 0, or by utilising the shower itself to generate the correct 
probabilities. This latter option is commonly referred to as the pseudoshower approach [HIH]- 
A complication arises if the phase-space separation criterion Q%j k is different from the parton-shower evolu- 
tion variable t. This can imply the possibility of a shower emission Q < Q cu t being allowed "between" two 
branchings at Q > Q cu t m the parton-shower history of the matrix element. In such cases, in order not to 
spoil the logarithmic accuracy of the parton shower, the existing branchings need to be embedded into the 
subsequent parton-shower evolution. This leads to a truncated shower algorithm [21 [T]. 



If the expectation value (O) in Eq. (2.16) is the total cross section, i.e. if O = 1, we can write 



(MEPS) 



J d$ s (M)B({a}) 

{/} 



A( ps )(to,^;{a}) 



E E 16tt 2 ,., 



dt 



dz 



unresolved 
2tt j i 



2tt 



JijAt,z,<l>)0{r~r{{a})) 



1 8tt a. 



Sij t 



(2.17) 



A( P W;{a}) 



£({a};t) 

x I 6 K? cut -Qtj,k(*,2,0)) +w{r~~ k (fi,t,z, <p;{a}))Q (Qij t k(t,z, 



resolved, PS domain 



resolved, ME domain 



where (cf. [TS]) 



w({d}) 



E E 

{m,n} Z^m,n 



S(b mn ,i({f})) B(b mn j({a})) 8^a s 



S({f}) 



K({3}) 



^mn,;({a}) 



(2.18) 



If w({a}) equals one, i.e. if the parton shower approximation is equal to the real-radiation matrix clement, 
the t- integral can be performed easily and the square bracket in Eq. ( 2.17| equals one. The total cross section 
is therefore identical to the leading order result, as a direct consequence of the unitarity constraint for the 
parton shower. The more common configuration will however be, that w({a}) ^ 1. In this case, the square 



bracket in Eq. (2.17) is different from one and the total cross section is not equal to the leading order result. 



The origin of the mismatch is a difference in the emission rate of the parton shower compared to the ratio 
R/B of matrix elements. While the former is exponentiated into the Sudakov form factor, the latter appears 
in the differential real-radiation probability only. We refer to this effect as an "emission-rate difference" in 
the following. 



To investigate its consequences on arbitrary observables, Eq. (2.16) can be expanded in powers ofa s , resulting 
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(O) 



(MEPS) 



= W d$ B m) B({3» 0({a}) 
{/} 



EE E 

{/} M.fc} 

87r a s 1 

EE E 

{/} ivM fi=9>9 

87r a, 1 



1 r d* B m) f dt 

J to 



dz 



2jv 



2tt 



B({a})/Cy,fe(t,z,< 



A^,k(/i>M,<Ma});*) 



C({a};t) 



2tt 



27T 



Jij,k{t, z, 4>) 

0(r^~ k ({a}))-0({a}) 

Jij,k(t, z, 4>) 



£( r ij,fc(/i>M>0; {#});*) 



, B({a}) JC^ft, s, 0) " v — f - 0(r r ,i({a})) 



(2.19) 



£({a};t) 

x 6 \Qij,k(t,z,4>) - Q cut ) w(r~ ^(/i,t,2,0;{a})) - 1 
+ O(c^) ■ 

While the first and second terms represent the usual parton-shower prediction, the third term encodes the 
emission-rate differences effected by matrix element corrections in the region of well separated partons. In 
the ME+PS approach they enter at 0(a s ). 

Emission-rate differences seem to be an undesirable side-effect of the ME+PS method at first. However, 



given Eq. (2.191, they can serve as an indicator for the relevance of higher-order real-emission corrections. 



We will elaborate on this fact in some more detail in Sec. BJ 



3 Merging Powheg and ME+PS - The MEnloPS approach 



In this section, the two master equations for the POWHEG (Eq. (2.11 )) and ME+PS (Eq. (2.16)) approaches 
are combined into one single expression, defining the MENLOPS approach. The aim of this combination 
algorithm is to simultaneously have NLO accuracy in the cross section, leading logarithmic accuracy as 
implemented in the parton shower and hard higher-order emissions corrected using tree-level matrix elements. 



Our method of choice is to simply replace the unresolved and the PS resolved part in Eq. (2.16) with the 



respective POWHEG expression. This essentially amounts to the replacement of the parton-shower no-emission 
probability with the corresponding POWHEG result, A^ PS ^ — » A^ ME ) and a substitution of the leading-order 
weight B by B, like in the POWHEG method itself. 

The ME part of the cross section is then generated separately, starting from real-emission matrix elements, 
as described in Sec. |2.2| This immediately implies that it will not automatically benefit from a POWHEG 
implementation regarding the local if-factor B/B, and it is therefore necessary to supply this A"-factor 
explicitly. There is no a-priori definition of a Born- level parton configuration in this context, because the 
ME event is defined in terms of a real-emission configuration. One rather has to identify a branching history 
{a} — > foij,fc({a}) such that B/B can be computed depending on b l3 ^{{a})- The definition is achieved by 
clustering the real-emission configuration using an algorithm which is similar to a sequential recombination 
jet scheme and which determines the node to be clustered according to the related branching probability in 
the parton shower. For more details on this technique we refer the reader to Q]. 
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Implementing these ideas, the master formula for the first emission in MENLOPS is obtained as 



(O) 



(MENLOPS) 



J d<M{p})B({a}) 

{/} 



A<> ME \t ^ 2 ;{a}) 0({a}) 



unresolved 



E E 



{ij,k} fi=q,g 



16tt 2 



dt 



dz 



Jij, k {t,z,<f>)0{r~^{{(i})) 



1 S(r~ k({f})) %, fc (r ra ({a})) 



s » s({f}) ma}) u 



( A( ME )(i, / i 2 ;{a})e(Q cut -Q lJ , fe 



(3.1) 



resolved, PS domain 



A( PS )(i, Ai 2 ;{a})e(Q u - fe -Q c 



olvcd, ME domain 



Note that the arguments of Qu k have been suppressed for ease of notation. The resolved-ME part of the 
expression in square brackets exhibits an additional factor 



A( PS )(f, M 2 ;{a}) _ i + y, y, 1 



A( ME )(t,/i 2 ;{a}) 
1 



16tt 2 



dt 



dz 



{rj,ic} fi=q,g 
S(r tj ,- k ({f})) R^(^({a})) Sna 



— J ij>k (t,z,<t>) 
o 2lT 



S({f}) B({5}) t 



JCij k(t, z, 4>) 



C{r~ ;k {{a})-t) 



£({a};t) 



(3.2) 



compared to the POWHEG master formula. This makes the emission-rate difference in the MENLOPS method 
explicit. However, the expectation value of O is still determined correct to 0(a a ), as can be seen by explicitly 
expanding Eq. (|3.1| in powers of a s : 



(O) 



(MENLOPS) 



E / d<MM) (B+V + l) ({a}) 0({a}) 
{/} 



EE E /d*B(w)d*«i 



{/}{*l,k}A=9>9" 



R^(r tl ,fc({«})) ©(*o- t) - Sij, fc (r~ fe({a})) 



0({a}) 



E E E / *MM) d$«j* Rti,*(r,j iS ({3})) 0(t - to) e (Q cut - Q ij>k ) 0{r~^{{a})) (3.3) 

{/}{Q,fc}/<=9'9 

E E E / d< MM) d*«J* Ry,fc(r- £ ({a})) 6(* - to) 6 - Q cut ) 0(r~ :k ({a})) 

{/}{y,fe}A=?'fl 



Comparing this with Eq. (2.12), we find that the 0(a s ) accuracy of the MENLOPS method is identical 
to what is obtained from POWHEG. The potential mismatch between exact higher-order tree-level matrix 
elements and their respective parton-shower approximation, leading to a difference between A^ ME ) and A^ PS \ 
contributes terms of 0(a 2 s ) or higher as long as (Qij,k — Qcut) enforces t > to- The precise value of Q cu t 
must therefore be chosen such that this constraint is satisfied. 



Apart from being of order a s , the term in the square brackets of Eq. (3.2 1 should be rather small in practice, 



as potential differences between matrix-element and parton-shower expressions merely lie in subleading 
logarithmic and power corrections. This corresponds to saying that the mismatch between POWHEG and 
MENLOPS results is at most of order a 2 log (^VQcut) if the parton shower has LL accuracy, and of order 



if it has NLL accuracy. 



Generating a second emission using the ME+PS method supplemented with the above metioned local K- 
factor of course introduces additional emission-rate differences, as described by Eq. (2.19). However, because 
such terms are of 0(a 2 ), they do not spoil the next-to-leading order accuracy of the method. 



8 



At this point we would like to stress that in their publication Hamilton and Nason arrived at the same 
ideas [3]. 



4 Results 

This section collects results obtained with an implementation of the algorithm described in the previous 
sections in the SHERPA framework. It aims at detailing the improved description of data collected in vari- 
ous collider experiments and at quantifying some of the systematic uncertainties inherent to the MENLOPS 
method, in particular those related to the merging of the multijet tree-level contributions. Note again, that 
the MENLOPS approach is designed to merge the next-to-leading order accurate description of a given core 
interaction (like for example e + e~ — > qq) through the POWHEG method with higher-order tree-level contri- 
butions (like e + e~ — > qqgg) described in the ME+PS approach. Since the total cross section is essentially 
defined by the POWHEG expression of the core process in question, uncertainties like those related to the 
choice of scales are encoded mostly there. They have been discussed in our parallel publication |19j . while 
uncertainties related to the ME+PS method were discussed for example in [Tl 111], 

However, a comparison with results of the ME+PS and POWHEG techniques alone is extremely useful to 
assess the quality of the approach and the improvements related to it. The precise setup of SHERPA for 
this comparison, including in particular a parton shower based on C at ani- Seymour subtraction terms [22] 
and an automated implementation of the Catani-Seymour subtraction method |23J in the matrix-element 
generator AMEGIC++ [23] was described in detail in our parallel publication [T!5]. Throughout the paper, 
we use the CTEQ6.6 parton distribution functions [25], (and, correspondingly, the MS subtraction scheme) 
with as(mz) = 0.118 and running at two-loop. If not stated otherwise, hadronisation is not accounted for. 
Multiple parton interactions are not included in the simulation. We use the Rivet program package |26) and 
the HZTool library [27] for analyses and comparison with data. 



Merging Systematics 

As pointed out in the previous section, the ME+PS approach violates the unitarity of the parton-shower 
simulation. This discrepancy is directly inherited by the MENLOPS method. The extent of this effect depends 



entirely on the quality of the parton-shower algorithm, as can be seen in Eq. (3.2): If the parton-shower 
approximation to the real-emission matrix element is good, the correction factor, Eq. (3.2 1, is close to one. 

We test the quality of our algorithms in the reactions e + e~ — » hadrons, deep-inelastic lepton-nucleon scatter- 
ing, Drell-Yan lepton-pair production, W- and Higgs-boson production and W + W~ -production by varying 
the phase-space separation cut, Q cu t, and the maximum number of partons, A max , which is simulated with 
matrix elements in the MENLOPS approach. This is in close correspondence to the sanity checks of the 
ME+PS method which have been presented in PQ. The respective results are summarised in Tabs. [l][6] and 
depicted in Fig. [T] It is interesting to note that differences in the total cross-section are smallest for the 
MENLOPS samples with A max = 1, and that they increase steadily for larger iV m ax. This indicates that 
the parton shower tends to underestimate the cross section of higher-order tree-level contributions. We 
observe two important effects, which allow us to judge the quality of the MENLOPS approach with respect 
to NLO accuracy: Firstly, for A^ ma x = 1 the emission-rate differences never exceed the size of the NLO 
corrections. Secondly, for any given value of Q C ut and iVjet, the relative difference between the cross sections 
from MENLOPS and POWHEG is always smaller than the one between the cross sections from ME+PS and 
LO+PS. This is best seen in Fig. [T] and it gives some confidence that the MENLOPS technique can help to 
improve perturbative QCD predictions from parton-shower Monte Carlo. 

The above analysis can be seen from a different perspective as follows: Usually, the biggest intrinsic uncer- 
tainty of the ME+PS approach stems from the freedom to choose the phase-space separation cut, Q cu t, as 
explained and exemplified in a number of processes in pQ. Since the MENLOPS method relies on identical 
ideas to separate the real-emission phase space, it naturally inherits this source of uncertainty. Deviations 
of MENLOPS results from results with different values of Q cu t are to be expected. However, their small size 
in a reasonable range of Qcut is a sign of the algorithm working well. The following rule of thumb can be 
applied: If the value of Q cut is chosen too large, too much extra emission phase space is left to the POWHEG 
simulation, typically leading to an underestimation of jet rates, since POWHEG only simulates the first emis- 
sion through matrix elements. If, on the other hand, this value is too small, too much phase space is filled 
by matrix elements with large final-state multiplicity, which may lead to noticeable emission-rate differences. 
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logio y cut 


max 





1 


2 


3 


LO+PS 


POWHEG 


ME+PS 


MEnloPS 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


-1.25 


28.37(1) 


29.43(1) 


28.09(2) 


29.60(2) 


28.11(2) 


29.63(2) 


28.09(2) 


29.63(2) 


-1.75 


27.46(2) 


29.47(3) 


27.56(3) 


29.43(3) 


27.48(3) 


29.43(3) 


-2.25 


27.11(3) 


29.38(3) 


26.90(3) 


29.18(4) 


26.93(4) 


29.17(4) 



Table 1: Dependence of the inclusive e + e — > jets cross section in nb on the number of extra jets generated 
in both the ME+PS and the MENLOPS method. For N max = these reduce to the LO+PS and 
POWHEG methods, respectively. 



Qcut 


N 

1 y max 





1 


2 


3 


LO+PS 


POWHEG 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


3 GeV 


365.7(4) 


322.4(6) 


380.5(6) 


333.5(9) 


388.4(6) 


339.3(1.0) 


396.6(9) 


344.4(1.2) 


5 GeV 


380.0(6) 


330.6(9) 


387.0(6) 


339.2(1.1) 


393.2(8) 


342.8(1.2) 


9 GeV 


381.8(7) 


330.4(1.0) 


388.1(8) 


337.9(1.2) 


394.2(8) 


344.0(1.2) 



Table 2: Dependence of the DIS production cross section in nb on the number of extra jets generated in 
both the ME+PS and the MENLOPS method. For iV max = these reduce to the LO+PS and 
POWHEG methods, respectively. The precise definition of Qcut is given in [IT] . 



Qcut 


max 





1 


2 


3 


LO+PS 


POWHEG 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


15 GeV 


1993(1) 


2423(1) 


2114(6) 


2549(9) 


2169(6) 


2587(9) 


2159(7) 


2599(10) 


20 GeV 


2103(4) 


2516(6) 


2137(5) 


2548(9) 


2135(5) 


2548(9) 


40 GeV 


2092(3) 


2477(9) 


2104(3) 


2485(7) 


2101(3) 


2482(7) 



Table 3: Dependence of the W production cross section in pb on the number of extra jets generated in both 
the ME+PS and the MENLOPS method. For V max = these reduce to the LO+PS and POWHEG 
methods, respectively. 



Qcut 


max 





1 


2 


3 


LO+PS 


POWHEG 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


15 GeV 


394.7(1) 


477.8(1) 


417.7(8) 


489.3(8) 


425.9(9) 


502.6(1.3) 


429.4(1.1) 


503.6(1.3) 


20 GeV 


416.8(7) 


487.2(9) 


424.3(8) 


496.4(1.0) 


423.6(8) 


496.1(1.0) 


40 GeV 


417.1(4) 


486.6(6) 


419.8(5) 


489.1(6) 


420.6(5) 


489.1(6) 



Table 4: Dependence of the Z production cross section in pb on the number of extra jets generated in both 
the ME+PS and the MENLOPS method. For N max = these reduce to the LO+PS and POWHEG 
methods, respectively. 



Qcut 


max 





1 


2 


3 


LO+PS 


POWHEG 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


15 GeV 


1.063(1) 


2.425(1) 


1.217(4) 


2.698(11) 


1.297(6) 


2.846(15) 


1.353(11) 


2.937(27) 


20 GeV 


1.195(3) 


2.627(9) 


1.255(5) 


2.762(13) 


1.270(7) 


2.769(19) 


40 GeV 


1.177(2) 


2.488(6) 


1.215(3) 


2.571(10) 


1.212(4) 


2.598(13) 



Table 5: Dependence of the Higgs production cross section in pb in gluon fusion on the number of extra 
jets generated in both the ME+PS and the MENLOPS method. For AT max = these reduce to the 
LO+PS and POWHEG methods, respectively. 
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QcUt 


max 





1 


2 


3 


LO+PS 


POWHEG 


ME+PS 


MEnloPS 


ME+PS 


MENLOPS 


ME+PS 


MENLOPS 


20 GcV 


0.949(1) 


1.333(1) 


1.305(4) 


1.369(5) 


1.391(6) 


1.469(7) 


1.405(7) 


1.470(9) 


40 GcV 


1.277(3) 


1.379(5) 


1.361(4) 


1.459(6) 


1.378(6) 


1.465(6) 


80 GcV 


1.227(2) 


1.405(5) 


1.285(3) 


1.451(5) 


1.298(5) 


1.460(5) 



Table 6: Dependence of the W + W -production cross section in pb on the number of extra jets generated 
in both the ME+PS and the MENLOPS method. For N max = these reduce to the LO+PS and 
POWHEG methods, respectively. 



e+e~ -> jets at y/s = 91.25 GeV 



+ 1.04 

v 

> 1.02 



0.98 
0.96 
0.94 



l0 § g#=-i-25/-i-75/-2-25 (ME+PS) 
logg|h-i.25/-l.75/-2.25 (MENLOPS) 



pp -> W+jets at y/s = 1.8 TeV 



1.14 
1.12 

1.1 
1.08 
1.06 
1.04 
1.02 

1 .0 
0.98 



— Q cut =i5/20/40GeV (ME+PS) 
o- Qcut=i5/2o/4oGeV (MENLOPS) 



pp — > h+jets at y's = 14 TeV 



t 1-35 + 

eg z 

^ !-3 + 
b 

1.25 =- 

1.2 + 

1-15 + 

1.1 + 

1.05 + 

1.0 = — 

0.95 e_ 



Qcut=i5/2o/40 GeV (ME+PS) 
■o- Qcut=i5/2o/40 GeV (MENLOPS) 



e+p -> e++jets at = 300 GeV 



1.14 - 

1.12 — 

1.1 — 

1.08 + 

1.06 + 

1.04 + 

1.02 — 

1.0 - — 

0.98 — 



— Qcut=3/5/9 GeV (ME+PS) 
-°- Qcut=3/5/9 GeV (MENLOPS) 



=9= 



pp -> Z+jets at = 1-96 TeV 



1.14 - 

1.12 — 

1.1 — 

1.08 + 

1.06 + 

1.04 + 

1.02 — 

1.0 - — 

0.98 — 



Qcut=i5/20/40 GeV (ME+PS) 
-°- Qcut=i5/2o/ 4 o GeV (MENLOPS) 



pp -> W+W~+ jets at ,/s = 14 TeV 



t 1-7 

tT 1.6 

b 

1-5 
1.4 

1-3 
1.2 
1.1 
1.0 
0.9 







: -•- Qcut= 


=20/40/80 GeV (ME+PS) 


^ Qcut= 


=20/40/80 GeV (MENLOPS) 





















Figure 1: Visualisation of the emission-rate differences induced by the ME+PS and 
MENLOPS master formulas, Eqs. (3.1| and (2.16), respectively. 
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Figure 2: Differential jet rates d nn+ i for three different merging cuts, Q cu t, in Drell-Yan 
lepton-pair production at the Tevatron at y/s — 1.96 TeV. 
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The value of Q cut should therefore lie well between the parton-shower cutoff and the factorisation scale of 
the core process, with some margin on either side of this interval. 

We exemplify the stability of our MENLOPS implementation with respect to variations of <5 cut in Fig. [2] 
Due to their similarity to Qij,k, the differential jet rates shown there are extremely sensitive to the details of 
the radiation pattern and thus to the accuracy of the ME+PS implementation. They tend to expose even 
the slightest mismatch between PS and ME subsamples, which then shows up as a kink in the distribution. 
However, when varying Q cu t in a rather wide range, we observe no sizable discrepancies between the respective 
MENLOPS predictions, which is a very encouraging result regarding the quality of the algorithm and its 
implementation in SHERPA. 

4.1 e + e — > jets 

In this section we focus on electron-positron annihilation into hadrons at LEP energies (y/s =91.25 GeV). 
The core process of the simulation is therefore the reaction e + e~ — > qq. A full wealth of experimental data 
has been provided by the LEP experiments, which allows to assess the quality of the MENLOPS approach 
in this simplest realistic scenario. Although the improvements discussed in this paper concern only the 
perturbative QCD part of the Monte-Carlo simulation, our results account for hadronisation effects using 
the Lund model [301 E] to make them comparable to experimental data. Virtual matrix elements needed for 
the simulation were supplied by code provided by the BlackHat collaboration [31 j . 

Figure|3]highlights the improvement in the description of jet data. In the hard-emission region the MENLOPS 
results for the 2 — > 3-, the 3 — ?> 4- and the 4 — > 5-jet rate are generally closer to the data than the POWHEG 
ones, which hints at the success of the simulation. Deviations in the 5 — > 6-jet rate are most likely due 
to the fact that matrix elements for six-jet production are not included. Note that these distributions are 
normalised to the total cross section, such that no rate difference between the ME+PS and the MENLOPS 
samples can be observed. 

Figures [4] and [5] show examples of event-shape variables, which are all very well described in the hard-emission 
region by the MENLOPS simulation. Several distributions for jet angular correlations in 4-jet production, 
that have been important for the analysis of QCD and searches for physics beyond the Standard Model are 
investigated in Fig. [6j The good fit to those data proves that correlations amongst the final-state partons 
are correctly implemented by the higher-order matrix elements. 

4.2 Deep-inelastic lepton-nucleon scattering 

Deep-inelastic scattering (DIS) is one of the best understood processes in perturbative QCD. However, it has 
been an obstacle for a very long time to properly simulate hadronic final states in DIS using general-purpose 
Monte Carlo based on collinear factorisation. Only recently, a consistent approach was presented [IT], that 
allows to describe jet data throughout the experimentally accessible range of Q 2 , the negative virtuality 
of the exchanged virtual 7*/Z-boson. It is absolutely mandatory for this method that a large number of 
final-state partons can be described by hard matrix elements in order to lift the severe restrictions on the 
real-emission phase space of the parton shower, which are imposed by the factorisation theorem. 

In all our simulations we use the core process e + q — > e + q. We present results for two analyses. The first 
is the measurement of inclusive jet production in [32] . which covers different ranges of jet-pseudorapidity in 
the laboratory frame, r]i a t, in the low-Q 2 domain 5 < Q 2 < 100 GeV 2 . Jets are defined using the inclusive 
fey-algorithm 34J and are constrained to Et,b > 5 GeV and the pseudorapidity range —1 < < 2.8, where 
Et,b is the jet transverse energy in the Breit frame. The second analysis corresponds to the measurement 
of dijet production in 55], which covered a wider range of Q 2 and produced many doubly differential jet 
spectra. The acceptance region is 5 < Q 2 < 15000 GeV 2 and —1 < i]i a b < 2.5. Jet transverse energies are 
subject to the cuts £t,bi,2 > 5 GeV and Et,bi + Et,b2 > 17 GeV. The latter requirement is introduced 
to avoid Et,bi ~ £r,B2, which is the region of the phase space where next-to-leading order corrections are 
unstable due to implicit restrictions on soft emissions [55] , 

As outlined in [TTJ, a crucial observable is given by the inclusive jet cross section, differential with respect 
to E\ b /Q 2 . For E\ b /Q 2 > 1 it probes a part of the phase space where leading order Monte-Carlo models 
without the inclusion of low- a; effects are bound to fail in their description of jet spectra. Another very good 
observable to validate the proper Monte-Carlo simulation is the dijet cross section as a function of Q 2 . While 
still a relatively inclusive quantity, it is an important indicator for the correct simultaneous implementation 
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Durham jet resolution 3 — > 2 (Ecms = 91.2 GeV) 



Durham jet resolution 4 — > 3 (Ecms = 91.2 GeV) 





Durham jet resolution 5—^4 (Ecms = 91.2 GeV) 



Durham jet resolution 6—^5 (Ecms — 91-2 GeV) 





-In (1/45) 



Figure 3: Durham d n n +\ jet resolutions at LEP compared to data taken by the ALEPH experiment [28] 



Total jet broadening (Ecms = 91.2 GeV) 



Jet mass difference (Ecms = 91.2 GeV) 





Figure 4: Total jet broadening and jet mass difference at LEP compared to data taken by the ALEPH 
experiment [28] . 
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C-Parameter (£ CM S = 91-2 GeV) 



Thrust (Ecms = 91-2 GeV) 




ALEPH data 
MENLOPS (5-jet) 
ME+PS (5-jet) 
POWHEG 





ALEPH data 
MENLOPS (5-jet) 
ME+PS (5-jet) 
POWHEG 




0.6 0.65 0.7 0.75 0.8 0.85 o.g 0.95 



Figure 5: C parameter and thrust distribution at LEP compared to data taken by the ALEPH experi- 
ment 1281. 
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POWHEG 
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Figure 6: Angles between the leading (in energy) four jets defined using the Durham algorithm with y cut 
0.008. Results at the parton level are compared to data from the OPAL experiment [29 . 



15 




Figure 7: Left: The inclusive jet cross section as a function of Ej. B /Q 2 in bins of r]i a i,, compared to data 
from the HI collaboration |32j . E%, B is the jet transverse energy in the Breit frame, while rji a b 
denotes the jet rapidity in the laboratory frame. Right: The dijet cross section as a function of 
Q 2 in bins of Et,x + Et.2, compared to data from the HI collaboration |33j . 
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Figure 8: The dijet cross section as a function of the Bj0rken variable xb, compared to data from the HI 
collaboration l33l . 
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Figure 9: The dijet cross section as a function of the dijet mass rrijj, compared to data from the HI 
collaboration 
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of inclusive DIS and the additional production of hard QCD radiation. The high quality of the MENLOPS 
prediction for the two above observables is confirmed in Fig. [7| Discrepancies in the description of the 
E\ B /Q 2 -spectrum in the forward region can be attributed to the fact that the simulation is limited to three 
additional partons in the hard matrix elements. This restriction is imposed by the usage of the matrix- 
element generator AMEGIC++ [23]. Figures [8] and [9] exemplify again that the MENLOPS simulation correctly 
predicts multijet differential distributions in all regions of the phase space, while the POWHEG approach fails 
in the low-Q 2 domain. 
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4.3 Drell-Yan lepton-pair production 



Results for lepton-pair production through the Drell-Yan process are compared to data from the Tevatron 
at y/s = 1.96 TeV in Figs. 10p3 using the core process qq — > ££, where I = e, /j,. The invariant mass of the 



lepton pair was restricted to be within 66 < m«/GeV < 116 in the simulation. The MENLOPS and ME+PS 
samples use tree-level matrix elements up to Z + 3 jets with a merging cut of Q cu t = 20 GeV. Virtual matrix 
elements are provided by BlackHat (31]. The Z — > £ + £~ decay is corrected for QED next-to- leading order 
and soft-resummation effects using the Yennie-Frautschi-Suura (YFS) approach [36] . 

The Tevatron experiments provide a wealth of measurements sensitive to QCD corrections in Drell-Yan 
production. Fig. 10 shows the transverse momentum distribution of the lepton pair in two different analyses 
from the D0 experiment. The left hand plot displays a very recent analysis using the Z —> fi/i channel [37] to 
measure the Z-p± distribution normalised to the inclusive cross section. It requires muons with p± > 15 GeV 
in a mass window of 65 < m w /GeV < 115 and with |?7| < 1.7. The muon signal is corrected to the particle 
level including photons clustered in a cone of radius R — 0.2 around each lepton. The plot on the right hand 
side stems from an analysis in the electron channel 38J which uses Monte-Carlo models to correct the leptons 
for all acceptances including the pseudorapidity range and minimal transverse momentum. Here we display 
the peak region of the transverse momentum of forward Z bosons with \yz\ > 2. The agreement between 
all three approaches and the measurement is outstanding. In the bins at p± < 10 GeV non-perturbative 
effects like the intrinsic transverse momentum of partons in a proton might play a role. Related Monte-Carlo 
models in SttERPA could be tuned to reach an even better agreement. Still, the Monte-Carlo prediction lies 
within the experimental error band over the full range. 



Two more measurements from the D0 experiment are displayed in Fig. 11 The pseudorapidity of the Z 



boson [39] was measured in the electron channel requiring electrons with p± > 15 GeV in the mass window 
71 < TO co /GeV < 111. Again, all three Monte-Carlo approaches agree very well with the experimental data. 
The right hand plot shows the azimuthal correlation between the Z boson and the leading jet [3D]. This is 
a measurement in the muon channel with the same selection cuts as described above. The distribution has 
been normalised using the inclusive Z cross section and the comparison shows that the three approaches 
underestimate the total rate for Z+jet production with respect to inclusive Z production by approximately 
10%. This might hint at the need for NLO accuracy also in the Z+jet process. It is remarkable though that 
the inclusion of higher-order tree-level matrix elements significantly improves the shape of the distribution 
with respect to the POWHEG simulation. 

The observables presented so far are mainly sensitive to the correct description of the leading jet. For that 
reason even the POWHEG approach is well capable of providing sufficient accuracy in their prediction. We 
now proceed to observables sensitive to higher-order corrections. 

Figure 12 (left) shows the inclusive jet multiplicity [41] for jets constructed using the D0 improved legacy 
cone algorithm [42] with a cone radius of R = 0.5 and p± > 20 GeV. Jets were required to lie in |?7| < 2.5 and 
to be separated from the leptons by AR(£,jet) > 0.4. While POWHEG agrees with the data for the 7Vj Ct = 1 
bin it fails to predict the rate of events with more than one jet. The MENLOPS and ME+PS predictions 
impressively demonstrate the effect of higher-order corrections provided by tree-level matrix elements up to 
the third jet. They agree with the measurement within the error bands for 7Vj et = 2,3 and as expected fail to 
predict the correct four-jet rate because no matrix-element corrections have been applied at that multiplicity. 

Transverse momentum spectra of the three leading jets accompanying the Z boson were measured by D0 
in [33]. The distributions in Fig. 12 (right) and 13 are normalised to the inclusive cross section for Z 
production and the jets have been constructed using the same settings as in the multiplicity measurement. 
Both MENLOPS and ME+PS deliver a very good description of these spectra while POWHEG fails to describe 
the rate and shape for the second and third jet. 
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Figure 10: The transverse momentum of the reconstructed Z boson in Drell-Yan lepton-pair production 
at the Tevatron at \fi = 1.96 TcV. Experimental data stem from the D0 experiment [571 [55] 
and are described in the text. 
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Figure 11: Rapidity of the reconstructed Z boson [33] (left) and azimuthal separation of the boson and the 
leading jet [40] (right) in Drell-Yan lepton-pair production at the Tevatron at y/s = 1.96 TeV. 
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Inclusive jet multiplicity pT of ist jet (constrained electrons) 
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Figure 12: Inclusive jet multiplicity [41] (left) and transverse momentum of the leading jet [43 (right) in 
Z+jets events at the Tevatron at y/s = 1.96 TeV. 
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Figure 13: Transverse momentum of the second and third jet [43] in Z+jets events at the Tevatron at 
= 1-96 TeV. 
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4.4 W^+jets Production 



In this section we focus on the production of VF-bosons and their subsequent decay into an electron-neutrino 
pair at the Tevatron at y/s — 1.8 TeV. The core process of the Monte-Carlo simulation is therefore qq' — > tv. 
The separation criterion is set to (Jcut = 20 GeV and up to three extra jets are taken into account. The 
electron-neutrino pair is required to have an invariant mass of vri &l/ > 10 GeV. The W — » ev decay is corrected 
for QED next- to- leading order and soft-resummation effects using the YFS approach [36]. Virtual matrix 
elements are supplied by BlackHat |31j . 

The left panel of Fig. [14] displays the transverse momentum of the VF-boson as compared to data taken 
by the D0 collaboration [44] , while the right panel shows the exclusive jet multiplicity of fc^-clustered 
jets (D = 0.7) with at least 20 GeV. Although the event sample generated using the POWHEG technique 
only provides the best match to the central value of the data, all three event samples are well within the 
experimental uncertainties. On the other hand, already in the rate of single-jet events deviations between 
the POWHEG sample and both the MENLOPS and ME+PS samples are visible, with the latter two agreeing 
very well. Similarly, the POWHEG sample underestimates the amount of radiation into the central detector 



region, as exemplified in Fig. 15 The right panel of this figure shows that, since the POWHEG approach 
is capable of modelling the second hardest emission using the soft-collinear approximation of the parton 
shower only, its description of the angular separation of the the first two hardest jets is missing prominent 
features originating in the wide angle region. These features are of course present in the approaches having 
fixed-order matrix elements at their disposal. 



Figure 16 shows the differential jet rates doi, di2, ^23 and ^34 using the above /cj_-algorithm. While the first 
three of them, for the matrix-element merged samples, are described by matrix element to matrix element 
transitions, only the softer part of doi is described by such a transition for the POWHEG sample. The harder 
part of the doi receives corrections by matrix elements of higher jet multiplicity which are clustered into a 
single hard jet first. Of course, these corrections are missing in the POWHEG sample. Furthermore, d\2 is 
described by a matrix element to parton shower transition only in the POWHEG sample. Hence, it strongly 
underestimates the amount of hard wide-angle radiation. Similarly, both c?23 and are described by the 
parton shower only in the POWHEG sample, showing the same behaviour. It is worth noting that both 
the MENLOPS sample, implementing local if-factors, and the ME+PS, scaled by a global if-factor, agree 
within their respective statistical uncertainties over the whole range, indicating the well known fact of the 
approximate momentum independence of the virtual corrections to the leading order process. 
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Figure 14: Transverse momentum of the W, compared to data taken by the D0 collaboration [44], and 
the exclusive jet multiplicity in inclusive W production at the Tevatron at yS =1.8 TeV. 



Pseudorapidity of leading jet Separation between leading and second leading jet 




-4-2024 1 1.5 2 2.5 3 3.5 4 4.5 5 

(/(jet 1) AR(jet 1, jet 2) 



Figure 15: Pseudorapidity of the hardest jet and angular separation of the first two hardest jets in inclusive 
W production at the Tevatron at V S — 1.8 TeV. 
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Figure 16: Differential jet rates d nn+ i in W production at the Tevatron at y/s = 1.8 TeV. 
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4.5 Higgs boson production 



This section presents predictions for Higgs boson production via gluon fusion at nominal LHC energies of 
■s/s = 14 TeV. As NLO corrections to the core process gg — » h — » t + t~ are rather large, tremendous efforts 
have been made to perform fully differential calculations at NNLO [IS] and several predictions have been 
presented which merged such fixed-order results with resummation at next-to-next-to-leading logarithmic 
accuracy (46] . In this publication, we have no means for an improvement of the resummed calculation, 
instead we are restricted by the limitations of the parton-shower model. However, the systematic inclusion 
of higher-order tree-level matrix elements through the MENLOPS method can yield a significant improvement 
of existing NLO predictions, thus partially closing the gap between full NNLO predictions and Monte-Carlo 
results. It was shown, for example, in |47) that the predictions from ME+PS algorithms are often competitive 
to NNLO results if only the shape, not the normalisation, of observable distributions is concerned. 

In our simulations we set m# = 120 GeV and we include the decay h — > t + t~, however, the analysis 
focuses on the properties of QCD radiation associated with production of the Higgs boson. The invariant 
r-pair mass is restricted to 115 < m TT /GeV < 125 at the matrix-element level. Virtual matrix elements 
are implemented according to 4S|. The decay h — > t + t~ is corrected for QED soft-resummation and 
approximate next-to-leading order effects using the YFS approach 36J. 

Figure [17] shows the transverse momentum spectrum of the reconstructed Higgs boson. We observe that the 
POWHEG and MENLOPS samples are very consistent in the prediction of this rather inclusive observable. On 
the other hand, differences are observed in the results for individual jet transverse momentum spectra, cf. 



Fig. 18 They increase with jet multiplicity and with increasing transverse momentum, as can be expected, 
since the higher multiplicity jets are described by the uncorrected parton shower in the POWHEG method. 
Deviations are also found in the prediction of the dijet separation in r\ — space, which is shown in Fig. |19| 
However, it was previously found that the ME+PS result yields a prediction which is very similar to the 
NNLO result [IT] . This feature is naturally retained in the MENLOPS simulation. 



4.6 W + W~+jets Production 

In this section we present predictions for the production of the TU+f— > e + v e ] W~[-» [i~v^\ final state at 
nominal LHC energies of y/s = 14 TeV. The lepton-neutrino pairs are required to have an invariant mass of 
rag u > 10 GeV each. The W —> Iv decays are corrected for QED next-to-leading order and soft-resummation 
effects using the YFS approach 36 . Virtual matrix elements are supplied by MCFM [3U1[SD]. Again, this 
study focuses mainly on the properties of QCD radiation accompanying the diboson production process. Up 
to three additional jets at Q cut = 20 GeV are simulated in both, the MENLOPS and the ME+PS sample. 
It is known that high-multiplicity matrix elements in the ME+PS approach yield sizable effects on total 
event rates and shapes in this reaction [5T], a feature which is inherited by the MENLOPS method. Setting 
the phase-space separation criterion to a rather low value compared to the average partonic centre-of-mass 
energy will thus always lead to sizable emission-rate differences, which might be an indication of potentially 
large higher-order corrections. A similar effect was observed in a recent analysis of Z-boson pair production 
in association with a hard jet [52]. While the NLO corrections to this process are comparably small at 
Tevatron, they can be rather large at nominal LHC energies. Restricting the available final-state phase 
space by a jet veto, the corrections were again limited to smaller values, which makes the importance of 
the ZZ+2 jets final state explicit. As we include up to three additional jets in our simulation of W + W~ 
production, we observe similar effects. 

Figure [20] displays the invariant mass of the lepton pair and the scalar sum of transverse momenta of the jets, 
leptons and the missing transverse energy, Ht- While the former is described very well by the next-to- leading 
order calculation used in the POWHEG sample and receives only mild corrections from higher-order matrix 
elements, Ht receives sizable corrections at rather low values already. The reason for this is easily found in 
the sensitivity of Ht to any jet activity and thus to higher-order matrix element corrections of the parton 
shower. This can be seen in comparison to Fig. |23[ where a veto on additional jet activity was applied. We 
exemplify in Fig.[2l]that the ME+PS part of the MENLOPS simulation predicts significantly harder radiation 
than the POWHEG subsample. The corresponding corrections naturally amplify the deviations between the 
respective predictions of Ht- We show the impact of a jet veto on this distribution in the right panel of 

Fig. [231 

Figure [22] presents predictions for the azimuthal separation of the leptons and the two hardest jets. Again, 
the former receives only comparably small corrections, while higher-order matrix-element corrections have 
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Figure 17: Transverse momentum of the reconstructed Higgs boson in the gluon-fusion process at nominal 
LHC energies (14 TeV). 
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Figure 18: Transverse momentum of the first and second hardest jet in Higgs-boson production via gluon 
fusion at nominal LHC energies (14 TeV). 
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Figure 19: Pseudorapidity difference and angular separation of the first and second hardest jet in Higgs- 
boson production via gluon fusion at nominal LHC energies (14 TeV). 



2G 



Invariant mass of oppositely charged leptons HT 




o 50 100 150 200 250 300 10 2 10 3 

m e + [GeV] Hf [GeV] 



Figure 20: Invariant mass of the electron- muon pair (left) and Ht (right) in W + W production at nominal 
LHC energies (14 TeV). 
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Figure 21: Transverse momentum of the first and second hardest jet in W + W production at nominal 
LHC energies (14 TeV). 
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Figure 22: Azimuthal separation of the electron and the muon (left) and of the two hardest jets (right) in 
W + W~ production at nominal LHC energies (14 TeV). 
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Figure 23: Invariant mass of the electron- muon pair (left) and Ht (right) in W + W production at nominal 
LHC energies (14 TeV) after vetoing events with more than one jet with p T > 20 GeV. 
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Figure 24: Azimuthal separation of the electron and the muon (left) and transverse momentum of the 
hardest jet in W + W~ production at nominal LHC energies (14 TeV) after vetoing events with 
more than one jet with > 20 GeV. 
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large impact on the latter. This hints at the importance to include higher-order matrix elements in Monte- 
Carlo simulations of hadron-collider events if the hadronic centre-of-mass energy is large. The effect of a jet 
veto on the azimuthal separation of the leptons is shown in the left panel of Fig. [24] 
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5 Conclusions 



In this publication, a parallel development and independent implementation of the MENLOPS algorithm, 
first discussed in [3], has been presented. This new algorithm combines the so far most advanced methods to 
include higher-order corrections to a given core process: The POWHEG technique, which allows to produce 
inclusive samples for that process with next-to-leading order accuracy, and the ME+PS technique, which 
allows to generate inclusive samples with a leading-order cross section, but with the production of additional 
hard radiation corrected by higher-order tree-level matrix elements. 

Until the work of Hamilton and Nason |3| and the work presented here, these two approaches were considered 
orthogonal and thus used independent from each other, in the regime of their respective strengths and validity. 
With the recent efforts on combining them, the shortcomings of each method, i.e. the description of higher 
jet multiplicities in POWHEG and lack of the correct NLO cross section in ME+PS, have been expunged. 
We fully confirm the findings of Hamilton and Nason concerning both the formalism and the relative improve- 
ment in the simulation obtained through it. This is even more emphasised here, since the implementation of 
the MENLOPS method as presented in [3] seems to suffer from the choice of tools. As already indicated in the 
introduction, the omission of truncated showering in the program used to simulate the ME region may have 
caused a few of the uncertainties. We are convinced that in total, the superior quality of the ME+PS part 
of the simulation in SHERPA, including the truncated showering, are the only reason behind the improved 
simulation here - the formalism is identical in both publications. The drastically reduced uncertainties stress 
the great improvement by the MENLOPS method. 

Our results and the ones presented by Hamilton and Nason in fact show a significant improvement of 
many aspects of previous simulations in a variety of processes, including here e + e~~ annihilation to hadrons, 
hadronic final states in DIS, jets in association with single vector bosons and with vector boson pairs, and 
the production of Higgs bosons through gluon fusion. 

In the future, the description of many more processes with this combined NLO matching and multijet merging 
will become feasible. This is possible, because both the POWHEG and the ME+PS part of the implementation 
are fully automated in SHERPA. 

We would like to also point out that the methods developed so far will naturally serve as a starting point to 
promote the ME+PS idea to full NLO, in the sense that merging sequences of multijet matrix elements at 
NLO into one inclusive sample becomes feasible. A first attempt to achieve this from a somewhat different 
angle has been presented in [53] . 
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